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Abstract 

We investigate the energy behavior of the diffraction cone in the exclusive vector 
meson production in diffractive DIS within the /^-factorization approach. In our cal- 
culations, we make full use of fits to the unintegrated gluon structure functions ex- 
tracted recently from experimental data on F% p . Confirming early predictions, we ob- 
serve that shrinkage of the diffraction cone due to the slope of the Pomeron trajectory 
is significantly compensated by the anti-shrinkage behavior of the 7 — ► V transition. 
In order to match recent ZEUS data on the energy behavior of the diffraction slope, 
a'gfJJ/ip, exp.) = 0.115 ±0.018(stat.)~^QQ^(syst.) GeV -2 , we had to use an input value 
a'jp = 0.25 GeV -2 . We investigate the compensation effect in detail and give predictions 
for Q 2 -dependence of the rate of cone shrinkage for different vector mesons. 

1 Introduction 

The exclusive production of light and heavy vector mesons in diffractive DIS 

7 (*)p _^ Vp (1) 

is an ideal testing ground ^ El El of the color-dipole approach, [IJ El E], and of the 
fci-factorization approach, [3 El El- The two approaches are linked by the transverse 
Fourier-Bessel transform from the color dipole size space to the momentum space: the 
former is based on the color dipole cross section, while the key ingredient in the latter 
approach is the unintegrated gluon structure function T{x g , k 2 ) of the proton. Recently, 
the unintegrated gluon density was extracted from the experimental data on Fi p and 
presented in a form of simple, ready-to-use parametrizations for x g < 0.01 and through- 
out the entire phenomenologically relevant domain of gluon momenta k, ^Hj. These 
fits now allow one to put many of previous qualitative predictions of the color dipole 
or /^-factorization approaches to quantitative grounds. For example, they have already 
been successfully used in calculation of (virtual) photoproduction of dijets off proton or 
nuclei ^112]. In the vector meson production, these fits have helped understand the 
DGLAP factorization scale in the case of p meson production, confirming the earlier 
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ideas of [HJ • It seems now timely to reanalyze another issue that critically depends on 
the properties of the gluon content of the proton, namely, the shrinkage of the diffraction 
cone in the vector meson production. 



In simple Regge-type models with an approximately linear Pomeron trajectory, the 
(momentum transfer) t-dependence of the intercept 

da ( W'\ 4[aip( * ) ~ 1] , , 




forces growth of the slope b of the diffraction cone 



b(W) = b(W ) + 4c4> log 




) 
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the phenomenon called the shrinkage of the diffraction cone. In QCD, the Regge limit is 
described by the BFKL equation |13j . However, in the fixed coupling constant regime, 
the BFKL equation lacks any intrinsic dimensional scale and leads to the fixed value of 
the intercept. Inclusion of the running coupling constant splits the fixed cut on the j 
plane into an infinite sequence of moving Regge poles [12], each of them having its own 
non-zero a' . Within the color dipole BFKL approach ^S], the value of a'jp was obtained 
from the solution of color dipole BFKL equation for the slope ^0]; its value was shown 
to be nonzero ^7] and to depend both on the energy and the size of the color dipole. At 
very high energies, when the BFKL asymptotics would fully develop and the saturation 
effects would not yet come into play, the slope of the Pomeron trajectory would tend 
to a certain constant value independent of the dipole size. This asymptotic value is 
governed by the gluon screening radius and in [I^], under certain initial conditions, it 
was estimated to be a'-jp(asymp.) = 0.072 GeV -2 . At HERA, however, we still reside in 
the subasymptotic region, and the values of a'jp can be still W and Q 2 dependent. 

If we now return to the process (|T|). we find that the shrinkage of the cone due to 
the non-zero a'jp is only a part of the truth. Thanks to the factorization properties, 
the diffraction slope in the vector meson production can be decomposed into the beam, 
target and exchange contributions. A detailed analysis of each contribution performed 
in |18j showed that the 7 — > V transition (beam) contribution might possess substantial 
subasymptotic behavior. This energy behavior was expected to be of the anti- shrinkage 
type, i.e. the corresponding slope contribution was predicted to decrease with energy 
growth. This partially compensates the shrinkage of the exchange contribution to the 
slope, and, therefore, leads to the effective rate of the cone shrinkage observed in a vector 
meson V production a'^JV) smaller than a'jp. 

Quantifying this compensation requires the knowledge of the gluon density in the 
soft region 1 . At the time of publication ^H]> it was known rather poorly. Now, thanks to 
the fits obtained in ^0] > we have a much better understanding of how the gluon density 
behaves in the soft region. Therefore, it is possible now to reanalyze the anti-shrinkage 
properties of 7 —* V transition and check whether this compensation effect is important. 

This is done in the present work. On the basis of the fet-factorization approach, we 
calculate the differential cross section of the exclusive vector meson production in diffrac- 
tive DIS. We observe that the compensation effect is very important numerically. As we 
show below, the input parameter a'jp = 0.25 GeV -2 leads to the effective rate of the 

1 To be accurate, one should understand the words "soft gluon density" only as "the dipole-proton interaction 
in the soft region, parametrized in terms of gluon density" . 
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cone shrinkage in J/tp photoproduction production o^// ~ 0-12 GeV~ 2 . This reduction 
implies that much care should be taken when the vector meson production results are 
interpreted as a direct probe of the properties of the Pomeron. 

The structure of this paper is following. In Section 2 we briefly remind the basic for- 
mulae of the /^-factorization approach to the calculation of the vector meson production 
in diffractive DIS. By performing the small-f expansion of a generic helicity-conserving 
amplitude, we qualitatively study different contributions to the diffraction slope and dis- 
cuss the origin of the anti-shrinkage properties of the 7 — > V transition. We then proceed 
to numerical results, which are presented in Section 3. These results are discussed in 
Section 4, and Section 5 contains conclusions of this work. 

2 The k t - factorization predictions for the diffrac- 
tion slope: a qualitative analysis 

2.1 The basic amplitude 

The basic formulae for the vector meson production within the /^-factorization approach 
are well known (see details in |Sj). We denote the quark and gluon loop transverse 
momenta and the momentum transfer by k, K, and A, respectively (here, vector sign 
denotes transverse vectors). The fraction of the photon lightcone momentum carried 
by the quark is denoted by z, while the fractions of the proton light cone momentum 
carried by the two gluons are x\ and x 2 . With this notation, the imaginary part of the 
amplitude of reaction can be written in a compact form: 

C V y/Airatem f d 2 K 2 -> -* \ [ dzd ^ 1*1 T\ Tl\ \ \ f A\ 

ImA = s — 2 J -^asiq )T (xi, x 2 , Kx, k 2 ) j ^ _ ip v {z,k) ■ I(Ay,A 7 ) ■ (4) 

where = k± ^A. The helicity-dependent integrands I(Ay,A 7 ) have form 
I S (L,L) 
I s (T, T) 

I S (L,T) 



I b (T,L) 

where 

1 1 1 1 
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^ r + k r — K r + A/2 r — A/2 

1 " ~(r + K) 2 + e 2 ~ (r-£) 2 + e 2 (f+A/2) 2 + e 2 (f - A/2) 2 +e 2 ' 

with r = (1 — 2z)A/2 and e 2 = z(l — z)Q 2 + m 2 . Finally, the strong coupling constant 

is taken at q 2 = max[e 2 + k 2 , k 2 ]. 

In the absence of A — R correlations, and for a very asymmetric gluon pair, the off- 
forward gluon structure function T (x±, x 2 , k\, k 2 ) that enters (jlj) can be approximately 
related to the forward gluon density T{x g , k) via 

T(xi, x 2 <C xi, ki, k 2 ) « F(x g , k) ■ exp ^-i& 3IP A 2 j . 

Here x g ~ 0.41xi; the coefficient 0.41 is just a convenient representation of the off- forward 
to forward gluon structure function relation found in |19j . Numerical parametrizations 
of the forward unintegrated gluon density T(x g , k) for any practical values of x g and 
k 2 can be found in ^U]. The slope 631P contains contributions from the proton impact 
factor as well as from the Pomeron exchange; experimentally, it can be accessed in the 
high-mass elastic diffraction. 

The vector meson wave function ipv( z , k) describes the projection of the qq pair onto 
the physical vector meson. It is normalized so that the forward value of the vector meson 
formfactor is unity, and the free parameters are chosen to reproduce the experimentally 
observed value of the vector meson electronic decay width r(V — > e + e~). In what 
concerns the shape of the radial wave function, we followed a pragmatic strategy. We 
took a simple Ansatz for the wave function, namely, the oscillator type wave function and 
performed all calculations with it. In order to control the level of uncertainty, introduced 
by the particular choice of the wave function, we redid the calculations with another wave 
function Ansatz, namely, the Coulomb wave function, and compared the results. Since 
these two wave functions represent the two extremes (very compact and very broad wave 
functions that still lead to the same value of the electronic decay width), the difference 
observed should give a reliable estimate of the uncertainty. This difference is typically 
given by factor of 1.5 for absolute values of the cross sections, while in the observables 
that involve ratios of the cross sections (such as slopes, intercepts, etc.) this uncertainty 
is reduced. Details can be found in [8} l20|. 

Note also that when deriving we treated the vector mesons as IS wave states 
and used the pure S-wave spinorial structure 5 M instead of 7^, [2T] . 



2.2 The contributions to the diffraction slope 

The main feature of the \t\ = A 2 -dependence of the cross sections is a pronounced 
forward diffraction cone, which can be, at very small A 2 , parametrized by a single slope 
parameter b (see Section 3.1 for a detailed discussion on the definition of the slope). Let 
us understand the various contributions to the slope by studying the small-t expansion of 
a typical amplitude. Due to factorization properties, one can approximately decompose 
the slope to the target, exchange and the projectile contribution, [18] : 

b = bp^p + 6ip + 6 7 »_>y . (6) 

In principle, the presence of the helicity-flip amplitudes represents yet another source of 
the t-dependence. Since the single helicity-flip amplitudes are proportional to y/\t\, one 
can introduce correction to the slope as 

da/dt non . fHp + da/dt f up = A{t) + B ■ \t\ -» A(t) ■ e~ b f u ^ . 
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However, since the helicity-flip amplitudes are small, this correction never exceeds 2% of 
the value of the slope. 

The first term in © is an intrinsically soft quantity and is related to the wave function 
of the proton. It is not calculable within pQCD, but its magnitude can be estimated 
from the proton charge radius. In our calculations, we introduced the following elastic 
formfactor to the amplitudes 

f(a 2 ) = 2 ; with A = 1 GeV > ( 7 ) 

1 + A 2 /A 2 " 



which leads to the proton impact factor contribution to the slope 

b p ^p = 4 GeV~ 2 . (8) 

The second contribution in © arises from the f-dependence of the Pomeron intercept. 
In our calculations, it was parametrized as 

bjp = a'jp log ; a'jp = 0.25 GeV" 2 , x = 3.4 • 10~ 4 . (9) 

This dependence was ascribed both to hard and soft components of the unintegrated 
gluon structure function, see ^U]; more sophisticated parametrizations can be put forth 
as well. This parametrization was obtained by the requirement that we describe well the 
recent ZEUS data on the shrinkage of J /if) photoproduction We underline that, due 
to the compensation effect to be explained in a minute, this a'jp = 0.25 GeV -2 reduced 
to «e// ~ 0.12-r- 0.13 GeV -2 in the differential cross section, in agreement with the 
experimental data. 

Finally, the third term in (JHJ) originates from the photon to vector meson transition 
amplitude. This contribution possesses both Q 2 and energy dependence, which can be 
understood as follows. 

We start with the amplitude L — * L and consider function $2 at large enough values 
of e 2 so that k 2 can be safely neglected. Expand it at small A 2 and average over all 
possible directions of A. The result reads 



2£ 2 

$9 



E 2 (K 2 +E 2 ) (k 2 +£ 2 ) 3 

To the leading log Q 2 , it 2 «£ 2 , and one has 



e 2 -R 2 )+A[z 2 + {l-z) 2 \R 2 + "' ^ 



4e 2 4e 4 I 



(10) 
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s 

After performing the relevant integration, one obtains that the amplitude L — > L is 
proportional to 

2 G(x O 2 ) 2 G(x 2 )-nr K2 ^ HX9 ^ 2) 
Ql Ql Ql 2 Ql V 

Here G(x g , Q 2 L ) is the conventional gluon density, rjL = ((1 — 2z) 2 )i and fi 2 is an appro- 

2 

priately defined soft scale. The relevant hard scale Q L , which is linked to the scanning 
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radius of the color dipole approach r 2 «-> 1/Qj,, was investigated in detail in |DJ. Thus, 
one obtains the following contributions to the slope 

Ql f 2 G{x„Q\) 

The first term in (|1 1|) is a perturbative contribution. At Q 2 = 0, it should be of the order 
of several GeV -2 , but it quickly falls off with the Q 2 growth. The second contribution 
is a predominantly soft quantity, its dependence on Q 2 is weak. 

The T —* T amplitude can be analyzed in a similar way, and the qualitative results of 
this analysis remain the same as for the longitudinal case. This is why we use the above 
qualitative results even when discussing the properties of the transverse amplitudes. 

2.3 The compensation effect: sources of the cone anti- 
shrikage 

The total slope of the diffraction cone is given by the sum of all three contributions, (jfifl. 
Since the target contribution to the slope is taken constant, the energy behavior of the 
slope originates only from &ip and b^^y. 

The exchange contribution to the diffraction slope, 6ip, logarithmically grows with 
energy and leads to the well-known shrinkage of the diffraction cone. In a simple Regge- 
type models, the properties of the Pomeron coupling to hadrons are assumed to be energy 
independent, and this contribution is the only source of the energy dependence of the 
diffraction cone. In a more elaborate models, such as the A^-factorization approach, the 
beam contribution to the diffraction slope, b-y^v, also possesses the energy dependence. 

On the basis of the above qualitative analysis, one can identify two dinstinct sources 
of this energy behavior. First, even within the leading logQ 2 , the second term in (|1 1|) 
depends on energy. Indeed, both J^"(x 9 ,/i 2 ) and G(x g ,Q L ) rise with energy, but the 
exponents of their rise are different, see [TQj . For Q L >1 GeV 2 , the integrated glue 

taken at q\ rises with energy faster than the differential glue at the soft scale /x 2 , and 
this contribution to the slope will decrease with energy rise. This source of the energy 
behavior of 6 7 ^y was also discussed in [TH] . 

The second source of the b^^y energy behavior appears, if we consider (fTUj) beyond 
the leading log Q 2 . In this case denominators will contain Q L + k 2 instead of just Q L . 
Due to specific properties of the unintegrated gluon density, the typical values of k 2 grow 
with energy even at fixed Q L . This is clear from Fig.^ taken from [HI], where we showed, 
in a single plot, how the unintegrated gluon density J-(x g ,ti 2 ) changes with energy (or 
1/xg) growth. One sees that the relative weight of the large k 2 region increases as we 
move from x g = I CP 2 to x g = 10~ 4 . This proves that Q 2 L + k 2 increases — and the 
contribution to the slope again decreases — with the energy rise. This effect mostly 
relies on specific, yet unavoidable, properties of the gluon density. 

Both sources of the energy behavior of the beam contribution to the slope are of 
the anti- shrinkage type. They lead to the partial compensation of the diffraction cone 
shrinkage, and eventually result in a' e ^j < a'jp. 

It is interesting to note that one can, in principle, disentangle these two sources of 
the compensation effect. What one needs is the study of the /^-factorization prediction 
for the ultrahigh energy behavior of oc' e ft- Note that the second contribution to the 
compensation effect works at full strength only at x g > I0~ 4 . At smaller x g , i.e. at higher 



6 




Figure 1: The ^-dependence of the unintegrated gluon density J-(x g , k 2 ) (fit D-GRV) for 
several values of x g . The solid and dashed curves correspond, respectively, to the full gluon 
density and to the soft contribution only. 

energies, the /{-shape of the unintegrated gluon density practically does not change, 
and the growth of average values of k 2 with energy stops. Thus, in this region, the 
compensation effect should be only due to the first mechanism. 



3 Numerical results 
3.1 Definitions of the slope 

Before comparing the results of the /^-factorization predictions of the diffraction slope 
with the experimental data, we would like to discuss the definition of the slope itself. 
The literal definition of the local slope of the diffraction cone as a logarithmic derivative 
of the differential cross section is 

Since the differential cross section flattens as \t\ increases, the value of the local slope 
will decrease with \t\ growth. In Fig. [2J we show the /ct-factorization calculation of the 
local slope for the J/ip photoproduction within the region \t\ < 1 GeV 2 . 
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Figure 2: The momentum transfer dependence of the local slope ()12|) for the photoproduction 
of J/ip meson. The solid and dashed curves correspond to total energies of the 7p collision 
equal to 50 GeV and 200 GeV, respectively. 

The values of the diffraction slope obtained in the experiment are often the results of 
exponential law fits to the measured differential cross section performed within a certain 
t interval. This procedure approximately corresponds to a finite- difference slope defined 
as 



where t\ and t<i define the region of the experimental fit. When comparing our predictions 
with such data, we will use precisely this definition of the slope. 

In literature, other definitions of the slope can be encountered, such as ^^\t=o or 
l/(|i|). If the differential cross section were a pure exponential law, all these definitions 
would lead to the same values. However, in a realistic situation, the offset among them 
can reach ~ 1 4- 2 GeV -2 . This should be kept in mind when comparing different results 
of the slope. 

3.2 The energy dependence of the slope parameter 

We start our analysis with the J ftp photoproduction. The local slope parameter, pre- 
dicted from the fcj-factorization approach, has already been shown in Fig. |2] within the 
region \t\ < 1 GeV 2 for two values of the W 7p . One observes a steady growth of the slope 
as the energy increases, which leads to the shrinkage of the diffraction cone. Motivated 
by the Regge-model considerations, this growth is usually described by the following law: 



The quantity a4//> w hich we will call the rate of the diffraction cone shrinkage, is often 
assumed to be equal to the slope of the Pomeron trajectory a'jp. However, as we argued 




(13) 




(14) 
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in the previous section, there are grounds to expect that these two will differ from each 
other. 

The rate of the shrinkage o4// can be analyzed quantitatively on a plot of b as a 
function of energy W. As can be seen from Fig.|2J the rate of the cone shrinkage is roughly 
independent of the value of \t\, therefore one can expect that different definition of the 
slope will still produce similar o4//- In Fig. we plotted the ^-factorization calculations 
of the energy dependence of the diffraction slope b(W) in the J/ip photoproduction. In 
order to have control on the level of uncertainty introduced by the wave function Ansatz, 
we calculated b(W) for both the oscillator and the Coulomb wave functions. 
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Figure 3: The energy dependence of the diffraction slope in J/ip photoproduction. Experi- 
mental data points are from ZEUS 

The same Figure contains also the recent experimental data from ZEUS j^2]. Keep- 
ing in mind the discussion of the previous subsection, we attempted in our calcula- 
tions to match the experimental procedure of the diffraction slope evaluation. The kf 
factorization values of the slopes were calculated according to (|13|) with \t\\ = 0.1 GeV 2 
and |t2| = 0.7 GeV 2 . As we mentioned above, it is precisely these J/ip photoproduction 
experimental data that we use to fix the two free parameters of the diffraction slope Q. 
Therefore the agreement between the curves and the data points is nothing else but just 
the double-check of the consistency of our calculations. 

As can be seen from this Figure, the fcf-factorization calculations do not predict b(W) 
to be strictly linear function of log(W), as at the low energy end of the plot the curves 
flatten out. This is a consequence of the fact that the compensation effect discussed 
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above can be well energy dependent (see also next subsection). However, at W > 70 
GeV, the approximate linearity is indeed observed. Note also that J ftp photoproduction 
at low energies, W < 25 GeV, corresponds to x g > 0.01. This is the very edge of the x g 
region the gluon density fits were devised for. We prefer to eliminate potential problems 
with applicability of the gluon density parametrizations and focus on region W > 50 
GeV. If we now evaluate the shrinkage of the cone between W = 50 GeV and W = 250 
GeV, we find the /^-factorization values 

, | 0.127 GeV" 2 for the oscillator wave function , ,„ _ x 

e H ~ I 0.121 GeV" 2 for the Coulomb wave function , 

which are, of course, in agreement with the experimentally measured value |22| 

a' e// (exp.) = 0.115 ± 0.018(atat.)lo:oii {ay at.) GeV" 2 . (16) 

The key observation here is that the o4// va l ue predicted by the fcf-factorization is 
sigfinicantly less than the input parameter a'jp = 0.25 GeV" 2 . 
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Figure 4: The t-dependence of the effective intercept a e ff(t) calculated from the energy growth 
of the differential cross section of the J /if) photoptoduction between points W = 50 GeV and 
W = 250 GeV. Experimental data points are from ZEUS, 22 . 

Another look at how the diffraction cone behaves with the energy growth is given by 
the energy rise of the differential cross section itself. Introducing the effective intercept 



10 



for the differential cross section 
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one can study how this intercept changes with t. The results of this study, together 
with ZEUS experimental data [22]) are shown in Fig. |1J Although the experimental 
data were available down to \t\ = 1.34 GeV 2 , we limited ourselves only to \t\ < 1 GeV 2 
region. Within this region, the ^-factorization calculations based on either of the two 
wave functions are in a qualitative agreement with the data. Besides, one sees that the 
shape of the vector meson wave function has rather minor effect on the diffraction cone 
shrinkage, and, in the subsequent plots we will give only results based on the oscillator 
wave function. 



3.3 The Q 2 and W behavior of the compensation effect 

It is known experimentally, and it follows from our analysis in Section 2 as well, that the 
values of the diffraction slope b are sensitive to the virtuality Q 2 and to the mass of the 
vector meson produced my. It is, therefore, interesting to check how the values of Q 2 
and my will affect the rate of the diffraction cone shrinkage a^. 
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Figure 5: The effective rate of the shrinkage of the diffractive cone a' e ^ for p meson and J/ip 
meson production as functions of Q 2 . The fcj-factorization results are calculated between W = 
50 GeV and W = 250 GeV, and for the oscillator wave function only. The photoproduction 
data points are from ZEUS, [221 1221 

In Fig. we show the /cj-factorization predictions of the Q 2 behavior of a' e ^{p) 
and a' e ff(J/i/j). We also show the two ZEUS data points available [231 [221, both of 
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which correspond to the photoproduction limit. What one sees on both plots is a slight 
decrease of a^// f rom ~ 0.13 GeV -2 down to ~ 0.10 GeV -2 , which takes place at typical 
Q 2 ~ niy. This rather weak sensitivity of a4// t° the values of Q 2 + my is in accordance 
with where, too, the value of a^// was found to depend mostly on energy but not 
on the scanning radius, which is related to the virtuality, or to be more precise, to the 
hard scale Q . 

As mentioned in the previous section, the ultrahigh energy behavior of «e// can hdp 
disentangle the two sources of the compensation effect. If jp collision energy grows to 
high enough values, such that x g < 10 -4 (which corresponds, for J/ip photoproduction, 
to W > 300 GeV) , then the second source of the anti-shrinkage behavior discussed in 
Sect. ESI should get suppressed. As a result, what one probes in this region is only the 
first contribution to the compensation effect. The overall compensation should, therefore, 
weaken. 

Numerical study shows that at W 3> 100 GeV the effective rate of the diffraction cone 
shrinkage a^yy significantly rises, in accordance with expectations. At W = 1000 GeV 
and W = 4000 GeV, a' e/ .rises up to 0.19 GeV" 2 and 0.22 GeV" 2 , respectively, which 
should be compared to (|15|) at HERA energies. This is, first, in agreement with our 
expectation that the compensation effect should decrease with energy growth. Second, 
this results can be treated as an evidence that, at HERA energies, both sources of the 
anti-shrinkage effect are of comparable importance. 

We would like to stress that this analysis is just a self-contained investigation of the 
energy behavior of the compensation effect, which is the property the fci-factorization 
approach. We do not attempt to predict the energy behavior of o4//- Such a prediction 
might be feasible only with a reliable information on the energy behavior of the input 
parameter a'jp. In our calculations, we used constant value of a'jp = 0.25 GeV -2 for all 
energies, which might be very inaccurate beyond the HERA energy range. 

4 Discussion 

4.1 Is this the Pomeron trajectory? 

The experimental results presented in Fig. |1] were explicitly interpreted in 22. as a 
measurement of the Pomeron trajectory. Here, we argue that this interpretation is, at 
least, not that straightforward. 

The most direct argument comes from the value of «e//- The numerical calcula- 
tions show that the rate of the diffraction cone shrinkage o^// predicted by the kf 
factorization approach is sigfinicantly less than the input parameter a'jp , which quantifies 
the i-dependence of the intercept of the underlying Pomeron exchange. Therefore, there 
exists a mechanism of compensation of the cone shrinkage. The qualitative analysis of 
Section 2 suggests that the energy decreasing contribution (fTT|) of the 7 — > V transition 
to the overall diffraction slope is in charge of this compensation effect. In Sect. |2~3*1 we 
pointed out two sources of this anti-shrinkage behavior, and the subsequent analysis gave 
evidences that both of them were equally important. 

In order to make sure that there is no other source for anti-shrinkage but the 7 — ► V 
transition, we made a double-check and switched off the A-dependence of integrands 
I s (Ay, A 7 ) in @. This made a^// j um P U P 

to ps 0.24 GeV" 2 , which is very close to the 
input value a'jp = 0.25 GeV -2 . This proves that it is precisely the 7 — > V transition 
vertex that causes such a strong reduction of the rate of the diffraction cone shrinkage 
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in vector meson production. 

Thus, we are led to a conclusion that, first, the experimentally measured values of 
a' e fj should not be interpreted directly as the slope of the Pomeron trajectory. Second, 
as we checked, o4// _a ip relation is rather robust within the /^-factorization scheme, 

a 'eff ~ a TP ~ O-l^ GeV~ 2 , 

and the experimental data shown in Fig. 0] might be in fact an evidence that the true 
value of the slope of the Pomeron trajectory, as measured in the J/ip photoproduction, 
is around 0.25 GeV -2 . 

It is interesting to note that this value is close to what is usually believed to be the 
soft 2 a', see [Ej- One can also quote a similar value measured in the elastic hadronic 
interaction at high energies |27l I28| . but, in its own turn, the single Pomeron exchange 
in these reactions receive sizable absorption corrections |2H] and, therefore, these results 
should not be interpreted straightforwardly as well. In any case, due to the strong 
compensation effect, we think that it would be premature to conclude that the measured 
value of ot-'eff in J/ip photoproduction is inconsistent with the soft Pomeron. 

4.2 Comments on other works 

We would like to mention that our approach to &' e ff is different from approaches of |25| 
and |26j . In the differential cross section of vector meson production was fitted to 
the experimental data with no attempt to disentangle the real physics that leads to this 
form of the cross section. An explicity non-linear Pomeron trajectory was introduced, 
which, in constrast to our predictions, leads to a strong t-dependence of a' . Although in 
this model a'(t = 0) = 0.25 GeV -2 , it becomes twice smaller at \t\ = 12m 2 « 0.24 GeV 2 , 
which allows the authors of [2>j to describe reasonably well the ZEUS data on shrinkage in 
J/ip photoproduction, [22]. It must be pointed out that this agreement is not surprising, 
since the fits to the cross section were derived precisely from the J/ip photoproduction 
data. By the construction of their model, the authors of directly relate the observed 
shrinkage in vector mesons photoproduction to the Pomeron properties, a4// = Q Jp- At 
the end, authors state that they "have reached a deeper understanding of the properties 
of dipole Pomeron", however, in the light of the present paper's results, such statement 
looks rather questionable. 

The authors of |26j worked in the color dipole formalism and, as can be expected 
from the early work ^S], they should have observed the compensation effect studied 
in the present paper. When discussing a', they indeed mention a similar effect, but 
they estimate it to be very inessential numerically, contrary to the claim of our work. 
The origin of this discrepancy lies in the fact that the authors of jSHj overlooked the 
main source of the energy dependence of the slope of 7 — > V transition and took only 
inessential part into account. Namely, in their phenomenological parametrization of the 
slope, the authors of |26j included only the first term in (jllj) (this is the (b 2 ) term in 
the notation of US])- The authors of 26_ neither took into account the second term in 
nor went beyond leading log Q 2 approximation by accounting for k 2 in this term, 
and, as a result, missed the sizable anti-shrinkage effect discussed here. Therefore, the 
conclusion a' e jj « a'jp of (22] was misleading. 

2 Ironically, several competing theoretical approaches agree on this, see also [241 1251 
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5 Conclusions 



In this work we investigated, at the quantitative level, the energy behavior of the diffrac- 
tion cone in the exclusive production of vector mesons. The work was conducted in the 
fc^-factorization scheme, closely related to the familiar color dipole formalism, and was 
based on recent fits to the unintegrated gluon density obtained in [10] . The fact that we 
do not devise models for the gluon content but instead heavily rely on the high-precision 
experimental data lends certain credence to the whole calculation. 

Confirming the results of early work ^H], we observed that the shrinkage of the 
diffraction cone, that could be expected from the Pomeron properties, is partially com- 
pensated by the anti-shrinkage behavior of the 7 — ► V transition vertex. We pointed 
out two sources of this compensation effects, both of equal importance, and observed 
that this reduction of the cone shrinkage is very significant numerically. In order to 
reproduce the experimentally observed value of a' e j^(J/tp) ~ 0.115 GeV~ 2 , we had to 
take ajp = 0.25 GeV -2 as an input parameter. This value turned out to be intriguingly 
close to what is usually believed to be a'jp of the soft Pomeron. 

This observation casts doubts on too straightforward interpretations of the experi- 
mental data on diffraction cone shrinkage as a direct measurement of the slope of the 
Pomeron trajectory. 

When studying the <5 2 -dependence of a^// within the ^-factorization approach, we 
observed some decrease of ot' e ff as we shifted from photoproduction limit to DIS. The 
exact numerical properties of this decrease depend on the particular choice of the defini- 
tion of the diffraction slope. This should not be forgotten when one compares the results 
of the theoretical predictions with the data. We also analyzed the energy behavior of the 
compensation effect and observed that it weakens with energy growth, which agree with 
expectations based on the qualitative analysis. Unfortunately, at this stage we cannot 
predict the energy behavior of the ct^// beyond the HERA energy range, since such a 
prediction requires understanding of the energy behavior of the input parameter ajp . 

In any case, the mere presence of the dependence of the compensation effect on kine- 
matics proves that the anti-shrinkage effect is not universal. This, in turn, supports the 
understanding that the whole underlying picture of the vector meson production (the 
"real" Pomeron) does not correspond to any simple singularity on j plane. 
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for his help at the early stage of the code development. I also wish to thank Prof. J.Speth 
for hospitality at the Institut fur Kernphysik, Forschungszentrum Jiilich. The work was 
supported by INTAS grants 00-00679 and 00-00366, and RFBR grant 02-02-17884, and 
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